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1. Introduction 

The currently prevailing view on the structure of the matter produced in nuclear 
collisions at the Relativistic Heavy-Ion Collider (RHIC) is anchored by two experimental 
observations [H El EJ HI |5] : (1) The emission of hadrons with a transverse momentum pt 
of several GeV/c or more is strongly suppressed of (jet-quenching), implying the presence 
of matter of with a very large color opacity, and (2) The anisotropic ("elliptic") flow 
in non-central collisions is near the ideal hydrodynamic limit, requiring an early onset 
of the period during which the expansion is governed by fluid dynamics (earlier than 1 
fm/c after the initial impact) as well as nearly ideal fluid properties with a viscosity-to- 
entropy density ratio r\js <C 1. The matter created at RHIC has been thus called the 
strongly interacting Quark-Gluon Plasma (sQGP) [6]. 

The origin of the jet-quenching phenomenon can be understood as follows [7J: 
during the early pre-equilibrium stage of the relativistic heavy-ion collision, scattering of 
partons causing the formation of deconfined quark-gluon matter often engenders large 
momentum transfers which leads to the formation of two back-to-back hard partons. 
These traverse the dense medium, losing energy and finally fragment into hadrons which 
are observed by the experiments. Within the framework of perturbative QCD, the 
process with largest energy loss of a fast parton is gluon radiation induced by collisions 
with the quasi-thermal medium. Elastic collisions add to the energy loss and are thought 
to be the dominant process for heavy quarks traversing the deconfined medium. 

Even though the phenomenon is being referred to as "jet-quenching", the 
overwhelming majority of computations of this effect have focused on the leading particle 
of the jet and do not take the evolution of the radiated quanta into account. A variety 
of schemes for quantitative calculations of the radiative energy loss have been developed 
[H [9j [TOj [HI [121 [13] . However, quantitative comparisons of theoretical model predictions 
incorporating a realistic hydrodynamic medium evolution with data from relativistic 
heavy ion collisions at RHIC [H] have revealed significant remaining ambiguities in the 
value of the extracted transport coefficients of the sQGP [15J. Different leading particle 
energy loss schemes may thus imply different values for these transport coefficients [7] . A 
systematic comparison of the assumptions underlying the various leading particle energy 
loss schemes and of their numerical implementations is currently under way [T6l [7]. 

Recently, with the advent of sophisticated experimental techniques for the 
reconstruction of full jets emitted from an ultra-relativistic heavy-ion collision [T71 ITS] , 
attention has shifted from leading particle energy-loss to the evolution of medium- 
modified jets. The study of the evolution of the entire jet in the medium is expected 
to lead to a better understanding of the dynamics of energy-deposition into medium 
and of the subsequent medium response, e.g. the possible formation of Mach-cones etc. 
[T9| l20l I2TJ . The current state of the art in for medium modified jets are Monte-Carlo 
generators [221 [231 [2U [25] , which calculate the medium-modified jet but do not take the 
medium response into account. The consistent treatment of the jet in medium as well as 
the medium response requires the application of transport theory, e.g. via Boltzmann 
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equation based calculation, such as is done in the Parton Cascade Model (PCM) [26J. 
Parton Cascades have already been applied to the time-evolution of ultra-relativistic 
heavy-ion collisions. However, for the most part the PCM calculations have focused 
on reaction dynamics [271 [28] . thermalization [29], electromagnetic probes [30, EI] and 
bulk properties of the medium [32], as well as to a far lesser extent on leading particle 
energy- loss [331 EH]- We would also like to note the recent progress in combining a 
Boltzmann equation based particle transport approach with pQCD cross sections for 
the scattering of high momentum particles with soft particle interactions mediated by 
a Yang-Mills field [35]. 

It is our long-term goal to advance the application of the PCM to the description 
medium modified jets and the respective response of the medium to a jet propagating 
through it and depositing energy in it. The achievement of this goal requires the 
validation of the PCM against analytic test cases which can be reliably calculated for a 
simplified medium, e.g. a pure gluon plasma or a quark-gluon plasma in thermal and 
chemical equilibrium. In the present manuscript we take a first step in this direction. We 
calculate elastic energy loss in an infinite, homogeneous medium at fixed temperature 
within the PCM approach and compare our results to analytic calculations of the same 
quantity. In addition, we calculate the rate of momentum broadening of a hard parton 
propagating through the medium and compare the results of our analysis to analytic 
expressions for the transport coefficient q. 

2. Quark Gluon Plasma and Parton Cascade Simulations 

The medium in our calculations is an ideal Quark-Gluon-Plasma ,i.e. a gas of u, d 
and s quarks and anti-quarks as well as gluons at fixed temperature T in full thermal 
and chemical equilibrium. In addition, we also conduct studies for a one-component 
gluon plasma in thermal and chemical equilibrium. For our transport calculation we 
define a box with periodic boundary conditions (to simulate infinite matter) and sample 
thermal quark and gluon distribution functions to generate an ensemble of particles at 
a given temperature and zero chemical potential. We then insert a hard probe, i.e. a 
high momentum parton, into the box and track its evolution through the medium. The 
medium particles may interact with the probe as well as with each other according to a 
Boltzmann Transport equation: 



where the collision term Ci is a nonlinear functional of the phase-space distribution 
function. Although the collision term, in principle, includes factors encoding the Bose- 
Einstein or Fermi-Dirac statistics of the partons, we neglect those effects here. 
The collision integrals have the form: 
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with 



and 



D(F k (x,p)) = U F k(^P) ~ U F k(^P) (3) 
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in, out 



Si is a statistical factor defined as Si = T7 K° nt \ with fQ n,out identical partons of 

species a in the initial or final state of the process, excluding the ith parton. 
The matrix elements \Ai\ 2 account for the following processes: 

99 ~^ 99 99 QQ 19 ->■ 99 

qq' —> qq' qq —> qq qg — > qg (5) 

(99 -> 99") (qq -»■ (99 ~> ?V) 

with g and q' denoting different quark flavors. The flavor changing processes in 
parenthesis are optional and can be disabled to study the effect of jet flavor conversion. 
The gluon radiation processes, e.g. gg — > ggg are not included in this study, but will be 
addressed in a forthcoming publication. The amplitudes for the above processes have 
been calculated in refs. [361 EZ] for massless quarks. The corresponding scattering cross 
sections are expressed in terms of spin- and colour-averaged amplitudes \M.\ 2 : 

For the transport calculation we also need the total cross section as a function of s which 
can be obtained from ([6]): 

/ (^] dQ 2 . (7) 

Since our medium is in full thermal and chemical equilibrium, we can use the 
effective thermal mass of a gluon and a quark in the system to regularize the cross 
sections |38|: 



H 2 D = 7ra s d p J 7^p|^/p(^ T )' ( 8 ) 

where d p is the degeneracy factor of a parton p and Ci is N c for gluons and (N 2 — l) / (2N C ) 
for quarks. Inserting the thermal distribution yields a Debye mass of a gluon of /i^ = gT 



in a thermal gluon system at temperature T and of [Md — y (2iV c + Nf)/6gT for quarks 
and gluons in a quark-gluon plasma. For example, the dominant elastic cross sections 
thus are: 

da 99 ^ 99 2 9 1 

~dqY~ ~ 4 (q 2 ± + /4) 2 ' (J) 

~W = 2m '(A^r- (10) 

do«~M ,4 1 
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where the first order of Feynman diagrams have been included. 

For our studies we use two distinct implementations of the Parton Cascade Model, 
the the Andong parton cascade code [39], and the VNI/BMS code [28]. Both codes 
have been modified to contain the same cross sections so that we can verify the outcome 
of our calculations through cross-comparisons between the two PCM implementations. 
For the sake of simplicity we keep the coupling constant fixed at a value of a s = 0.3. 

To obtain the temporal evolution of a high momentum parton propagating through 
a gluon or a quark-gluon plasma, a high energy gluon or quark of initial energy Eq is 
injected into or thermal (quark-)gluon-plasma box at the center of it's xy-plane with 
p = (0,0, p z = \fE~f). Following each individual scattering event involving our hard 
probe, we record the energies and momenta of the outgoing partons. The parton with 
the larger momentum in the final state is considered to be our continuing probe particle, 
due to the dominance of small angle scattering in the implemented cross sections. A 
simple estimate on the number of interactions the probe will undergo per unit time or 
length can be obtained via its mean free path in the medium: 

X s = — , (12) 
a T p 

where cy is the total thermal cross section and p the density of a medium. A gluon's 
mean free path in a gluon plasma is 1.1 fm at T = 300 MeV and 0.76 fm at T = 400 MeV, 
and in the quark-gluon plasma 0.83 fm and 0.59 fm, respectively. The quark mean free 
paths are about 9/4 times longer than those of gluon ones at same temperature. 

3. Results and Discussion 

In our work, we shall focus solely on the hard parton propagating through the medium 
and analyze its energy loss and momentum broadening as a function of time and distance 
traveled. We would like to point out that our calculation includes the full information 
on the medium-response as well, however, a detailed analysis on its characteristics it 
outside the scope of the present work and will be discussed in a forthcoming publication. 

Let us start by comparing the distribution of transverse momentum transfers 
experienced by gluons with Eq = 100 GeV and E = 500 GeV, respectively, while 
propagating through a gluon plasma at T = 300 MeV to the analytic expression given 
by the differential cross section (jUJ). 

Figure [1] shows that we find excellent agreement between our transport calculation 
and the analytic expression for Eq = 500 GeV. At the lower incident probe energy, 
Eq = 100 GeV we note a suppression in the high momentum transfer tail of the 
distribution which we can attribute to the effects of phase space and energy conservation, 
the gluon loses a substantial amount of energy while traversing a medium of 50 fm depth 
and thus for most of the time does not carry sufficient momentum to scatter with high 
momentum transfers while the analytic formula continuously assumes that the CM 
energy between the colliding particles is much larger than the Debye mass. 
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Figure 1. The Monte Carlo sampling of transverse momentum transfer at each 
scattering. The jet only is allowed to have a collision with medium particles. 



We now focus on the elastic energy loss of a high-momentum parton in medium. 
The left frame of Fig. [2] shows the energy as a function of time for a gluon with initial 
energy of 100 GeV (or 50 GeV, respectively), propagating through a quark-gluon plasma 
(QGP) or gluon plasma (GP) at a temperature of T=400 MeV. The right frame repeats 
the calculation for a light quark instead of a gluon. The calculation, which includes only 
elastic processes (i.e. the flavor exchange reaction channels have been disabled in order 
to unambiguously study the flavor dependence) clearly shows the anticipated linear 
decrease of jet energy with time. One should note, that analytic calculations usually 
study the jet energy as a function of distance traveled - in our case we substitute time 
for distance in order to have a quantity that is not affected by the periodic boundary 
conditions for the system in our calculations. We have verified the linear relationship 
between distance traveled and elapsed time with a slope near unity. 

The two frames clearly show the difference between a GP and a QGP in terms of 
energy- loss at the same temperature. This difference is due to the significantly higher 
parton density in a QGP compared with a GP for the same temperature, resulting in a 
larger number of scattering partners for the hard probe to interact with. The difference 
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Figure 2. Left: the gluon jet energy as a function of traveling time in the medium 
at T — AOOMeV: The travelling distance is proportional to the time. GP is a gluon 
plasma and QGP a quark-gluon plasma and the simulation does not include gg — ¥ ggg 
process. Right: the quark jet energy as a function of time: The penetration length is 
proportional to the time. Plasma temperature T = AOOMeV. 




Figure 3. Left: gluon energy as a function of time for different initial values of the 
energy in a gluon plasma. Right: the same for a fixed initial energy and various 
medium temperatures. The solid lines represent an analytical calculation (see text for 
details). 



is less pronounced for the quark probe than the gluon probe, due to the difference in 
their interaction cross sections. We also observe that the linear decrease of energy as 
a function of time tapers off in the long time limit as the probe energy approaches 
the thermal regime. While at the probe energies studied here (between 50 GeV and 
400 GeV) this occurs at times not relevant in the context of an ultra-relativistic heavy- 
ion collision, for smaller probe energies frequently seen at RHIC, this may be a significant 
effect. 

Fig. [3] shows the energy and temperature dependence of the elastic energy loss of a 
hard gluon in a gluon plasma. The calculations (symbols) are compared to an analytical 
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Figure 4. Left: Energy as a function of distance for a gluon propagating through a 
GP and a QGP at the same temperature and at equivalent entropy-densities. A scaling 
of the energy-loss with the entropy-density is observed. Right: distribution of probe 
energies after traversing the medium for 5, 10, 50, 70 and 100 fm at a temperature of 
300 MeV. 



calculation [301121]: 

dE f d?k ^ m . r , 9 . „. da , . 

-Tt=lwf F>{k '' T) r q ^-^ )v v, (13) 

where v = E — E' is the energy difference between before and after collision. Utilizing 
the characteristics of our medium and the cross sections in our calculation, this equation 
can be discretized to a form which lends itself to a comparison with our calculation: 

E p (z) = E p (0) - z^f° In 

Using an iterative procedure we can calculate E p (z) for the initial gluon energies and 
temperatures used in Fig. |3] and compare them to our calculation. The agreement 
between the analytical calculation and the PCM is remarkable and serves as validation of 
the PCM framework. The validation of the PCM calculation in a well-controlled infinite 
matter calculation in full equilibrium at fixed temperature is of significant importance, 
since the PCM can easily be used to study realistic dynamic systems far off equilibrium, 
e.g. an ultra-relativistic heavy-ion collision, which can only be poorly described in the 
framework of (semi-) analytic calculations. 

In Fig. [2] we observe a difference in the energy loss a parton suffers when 
propagating through a GP or a QGP. We understand this difference to be due to 
the different overall particle densities associated with a GP and a QGP at the same 
temperature. In order to validate this point, we initialized a GP and a QGP at 
an identical entropy density of s = 87.5 fm -3 , corresponding to a temperature of 
T = 457 MeV for a GP and T = 318 MeV for a QGP. The results of this calculation 
(also in comparison to a GP and a QGP at a temperature of T = 400 MeV) are shown in 



E p (z)T 



(14) 
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Fig. HJ for the same entropy density, both the GP and the QGP inflict a nearly identical 
amount of elastic energy loss to the hard probe. This scaling suggests that the entropy 
density, rather than the temperature, may be a robust quantity for the characterization 
of a thermal QCD medium by energy-loss measurements. We note that it is by no 
means clear whether the deconfmed medium in the early stages of a heavy-ion reaction 
more closely resembles a GP (or a liberated "glasma" [42] ) or a QGP in full chemical 
equilibrium [38J. Most likely the chemical composition of the deconfmed medium created 
in a ultra-relativistic heavy-ion collision changes significantly as a function of time, 
whereas its entropy density (after the initial equilibration) does not vary as much. 

The right frame of Fig. H] shows the distribution of probe energies for a gluon with 
initial energy of 50 GeV passing through a gluon plasma at temperature T = 300 MeV 
after 5, 10, 50, 70 and 100 fm, respectively. The gluon distribution at t = is a 
^-function at E = Eq. The Figure shows several interesting features. We find that 
elastic energy loss acts as a diffusion processes in momentum space, with the width of 
the distribution being a function of the traveling distance. Likewise the peak of the 
distribution shifts as a function of traveling distance. Even after traveling significant 
distances, about an order of magnitude larger than possible in a relativistic heavy-ion 
collision, the gluon has not thermalized with the medium, in which a medium particle 
has an average energy of about 0.9 GeV at T=300 MeV. 

Another important quantity for the characterization of the hot and dense QCD 
medium is the transport coefficient q, which is defined as [43] : 



with the squared momentum transfer q 2 = —t. Generally, q can be interpreted as the 
amount of squared transverse momentum per unit length a probe accumulates as it 
propagates through the QCD medium. This interpretation lends itself to a definition of 
q suitable for extraction from microscopic transport calculations: 



One should note that q in general should depend on the distance l z the probe has traveled 
through the medium, since the probe loses energy in the process, and thus the average 
momentum transfers in individual interactions of the probe with the medium may vary 
as a function of that distance. 

Figure [5] shows q for a gluon plasma as a function of distance traveled through 
the medium for several temperatures. The initial values of q depend strongly on the 
temperature, ranging from approximately 1.5 GeV 2 /fm at T = 300 MeV up to about 
9 GeV 2 /fm at 600 MeV temperature. For T = 400 MeV we find good agreement 
with recent calculations by [35] using a Boltzmann transport for hard binary collisions 
combined with soft interactions mediated by a collective Yang-Mills field and by [34] 
with the BAMPS parton cascade model in its binary scattering mode. At leading-log 
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Figure 5. Sum of squared transverse momentum transfers per unit length as a function 
of distance traveled for a 50 GeV gluon in a gluon plasma at different temperatures. 
This quantity corresponds to the transport coefficient q. 



order and in the weak coupling limit, using the same matrix element as in our PCM 
calculation, the following analytic expression has been derived for q [44J: 



q(A) 



a s T m 2 D In 



T 2 



m 



D , 



+ 4vr« 2 Ann 



(17) 



if the cut-off A > T. Comparing expression ffT7|) to expression (|14p and following the 
derivation of (fTTj) we find that the cut-off A 2 ~ E P T. Replacing A 2 by aE p (z)T (with a 
being a proportionality constant to be determined) in expression ( IT71) provides a good 
fit to our results for a choice of a = 0.1. 



4. Conclusions 

We have studied elastic energy loss of high energy parton in an infinite, homogeneous, 
thermal medium within the PCM approach and have compared our results to analytic 
calculations of the same quantity. In addition, we have calculated the rate of momentum 
broadening of a hard parton propagating through the medium and have compared the 
results of our analysis to an analytic expression for the transport coefficient q. We find 
good agreement between the PCM calculations and the analytic expressions (within the 
approximations used in both cases), giving us significant confidence that our transport 
approach provides a reliable description of a gas of quarks and gluons at temperatures 
above ~ 2Tc in the weak coupling limit. We expect that the results presented in this 
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manuscript can be used as a benchmark by other parton-based microscopic transport 
calculations. The validation of the PCM against the analytic test cases presented in this 
manuscript now allows us to advance the application of the PCM to the description of 
medium modified jets in relativistic heavy-ion collisions and the response of the medium 
to a hard parton propagating through it. 
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